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Abstract. The unexpected diversity of exoplanets includes a growing number of super-Earth 
, planets, i.e., exoplanets with masses smaller than 10 Earth masses. Unlike the larger exoplanets 
previously found, these smaller planets are more likely to have similar chemical and mineralogical 
composition to the Earth. We present a thermal evolution model for super-Earth planets to 
identify the sources and sinks of atmospheric carbon dioxide. The photosynthesis-sustaining 
habitable zone (pHZ) is determined by the limits of biological productivity on the planetary 
surface. We apply our model to calculate the habitability of the two super-Earths in the Gliese 
581 system. The super-Earth Gl 581c is clearly outside the pHZ, while Gl 581d is at the outer 
edge of the pHZ. Therefore, it could at least harbor some primitive forms of life. 
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5/5 ' 1. Introduction 

i Very recently, Udry et al. (2007) announced the detection of two super-Earth planets 

in the Gliese 581 system; namely, Gl 581c with a mass of 5.06 and a semi-major 
axis of 0.073 AU, and Gl 581d with 8.3 M ffi and 0.25 AU. Both mass estimates are 
minimum masses uncorrected for sini. The luminosity of Gl 581 can be estimated as 
L = 0.013 ± 0.002L© with a stellar temperature of T e = 3480 K and a stellar age of at 
, least 2 Gyr. The main question is whether any of the two super-Earths around Gl 581 
(T) ' can harbor life, i.e., that any of the planets lie within the habitable zone (HZ). Typically, 
stellar HZs are defined as regions around the central star, where the physical conditions 
■ are favorable for liquid water to exist at the planet's surface for a period of time long 
C > , enough for biological evolution to occur (Kasting et al. 1993). 

In the following, we adopt a definition of the HZ previously used by Franck et al. 
J> , (2000). Here habitability does not just depend on the parameters of the central star, but 
also on the properties of the planet. In particular, habitability is linked to the photosyn- 
thetic activity of the planet, which in turn depends on the planetary atmospheric CO2 
concentration, and is thus strongly influenced by the planetary dynamics. We call this 
definition the photosynthesis-sustaining habitable zone (pHZ). In principle, this leads to 
additional spatial and temporal limitations of habitability. 



2. Estimating the Habitability of a Super-Earth 

To assess the habitability of a super-Earth, i.e., a rocky planet smaller than 10 Earth 
masses (Valencia et al. 2006), an Earth-system model is applied to calculate the evolution 
of the temperature and atmospheric CO2 concentration. The numerical model couples 
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the stellar luminosity, the silicate-rock weathering rate, and the global energy balance to 
obtain estimates of the partial pressure of atmospheric carbon dioxide Pco 2 > the mean 
global surface temperature T sur f , and the biological productivity II as a function of time 
t. The main point is the persistent balance between the CO2 sink in the atmosphere- 
ocean system and the metamorphic (plate-tectonic) sources. This is expressed through 
the dimensionless quantities 

/wr(t) • f A (t) = /„(*), (2.1) 

where f WT {t) is the weathering rate, /a(£) is the continental area, and f ai (t) is the 
spreading rate normalized by their present values of Earth. The connection between 
the stellar parameters and the planetary climate can be formulated by using a radia- 
tion balance equation. The evolution of the surface temperature is derived directly from 
the stellar luminosity, the distance to the central star and the geophysical forcing ratio 
(GFR := fsr/fA)- For the investigation of a super-Earth under external forcing, we adopt 
a model planet with a prescribed continental area. The fraction of continental area to 
the total planetary surface is varied between 0.1 and 0.9. 

The thermal history and future of a super-Earth is determined by calculating the 
GFR values. Spreading rates can be derived from the mantle temperature. Assuming 
conservation of energy, the average mantle temperature T m is obtained as 

^pc(R 3 m - Rl)^ = -^R 2 m q m + T*E{t){Rl - R 3 C ), (2.2) 

where p is the density, c is the specific heat at constant pressure, q m is the heat flow 
from the mantle, E is the energy production rate by decay of radiogenic heat sources in 
the mantle per unit volume, and R m and R c are the outer and inner radii of the mantle, 
respectively. The photosynthesis-sustaining HZ (pHZ) is defined as the spatial domain 
encompassing all distances R from the central star where the biological productivity is 
greater than zero, i.e., 

pHZ := {R I n(P C o 3 (-R,t) j r Burf (iZ,*)) > 0}. (2.3) 

In our model, biological productivity is considered to be solely a function of the surface 
temperature and the CO2 partial pressure in the atmosphere. Our parameterization 
yields zero productivity for T surf < 0°C or T surf ^ 100°C or P C o 2 < !0 -5 bar (Franck 
et al. 2000). To calculate the spreading rates for a planet with several Earth masses, the 
planetary parameters have been adjusted following Valencia et al. (2006) as 



Rm ~ { M ffi 



(2.4) 



where R p is the planetary radius and M is the mass, with subscript © denoting Earth 
values. See von Bloh et al. (2007) for details and specified parameter values. 



3. Results and Discussion 

The habitable zone around Gl 581 for planets with five and eight Earth masses has 
been calculated for L = 0.013L Q . The results for Gl 581c and Gl 581d are shown in 
Fig. [T^,,b. The simulations have been carried out for a maximum CO2 pressure of 5 bar 
(light colors) and 10 bar (dark colors) neglecting the cooling effect of CO2 clouds. The 
super-Earth planet Gl 581c is found to be clearly outside the habitable zone. On the 
other hand, one might expect that life had a chance to originate on Gl 581d because it 
is near the outer edge of the pHZ, and for some combination of system parameters even 
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Figure 1. The pHz zone of Gl 581 for super-Earth Gl 581c (a) and Gl 581d (b) with a relative 
continental area varied from 0.1 to 0.9 and a fixed stellar luminosity of 0.013L©. The results are 
given as a function of planetary age. The light colors correspond to a maximum CO2 pressure 
of 5 bar, whereas the dark colors correspond to 10 bar. For comparison, the positions of Venus, 
Earth and Mars are shown scaled to the luminosity of Gl 581. 

inside the pHZ (von Bloh et al. 2007) . A planet of eight Earth masses has more volatiles 
than an Earth-size planet and can build up a sufficiently dense atmosphere to prevent it 
from freezing out due to tidal locking. Planets inside the habitable zone around M stars 
may be tidally locked, which however does not necessarily thwart habitability (Tarter 
et al. 2007). The modestly eccentric orbit of Gl 581d (e = 0.2 ± 0.1) further supports 
habitability, even if the maximum CO2 pressure is assumed to be as low as 5 bar. The 
appearance of complex life, however, is unlikely due to the rather adverse environmental 
conditions. To get an ultimate answer to the profound question of life on Gl 581d, we 
have to await the TPF/Darwin missions. They will allow for the first time the detection 
of biomarkers in the atmospheres of the two super-Earths around Gl 581. 
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